Quantum Memory Assisted Probing of Dynamical Spin Correlations 
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We propose a method to probe time dependent correlations of non trivial observables in many- 
body ultracold lattice gases. The scheme uses a quantum non-demolition matter-light interface, 
first, to map the observable of interest on the many body system into the light and, then, to store 
coherently such information into an external system acting as a quantum memory. Correlations of 
the observable at two (or more) instances of time are retrieved with a single final measurement that 
includes the readout of the quantum memory. Such method brings at reach the study of dynamics 
of many-body systems in and out of equilibrium by means of quantum memories in the field of 
quantum simulators. 



The quantum simulation of many-body physics with 
ultracold atoms in optical lattices requires the ability to 
prepare, manipulate, and probe the states [1]. Seminal 
experiments in this field range from the realization of a 
Bose-Einstein condensate with alkali atoms in the weak- 
coupling regime [2, 3], to the Mott-insulator-superfluid 
transition in the strongly interacting regime [4, 5]. Due 
to the high degree of control offered by ultracold lattice 
gases, simulation of quantum magnetism is becoming one 
of the most ambitious goals in the field [6-8]. Candidate 
systems are made up by optical lattices where the parti- 
cles (either bosons or fermions) arrange themselves in a 
deep Mott insulator, and the internal degrees of freedom 
of the atoms play the role of spins. Such internal states 
might be given from hyperfine structure in the case of 
alkalis [1] or by nuclear spins for alkaline-earth fermions 
[9]. At sufficiently low temperatures and entropies, the 
spin-spin interactions arising perturbatively from super- 
exchange processes (virtual tunneling between neighbor- 
ing sites) are predicted to give rise to quantum mag- 
netism phenomena as Neel ordering, SU(N) magnetism, 
spin Hall effects, and Stoner magnetism, see e.g. [10]. 

Quantum simulation with cold gases does not only 
aim at mimicking the phenomena encountered in con- 
densed matter physics, but also at exploring new fron- 
tiers in physics. In particular, ultracold atoms permits 
to explore non-equilibrium phenomena in closed systems 
[11]. Within this context, the study of thermalization 
in out-of-equilibrium situations after a quench, the infiu- 
ence of metastable states during the dynamics, the effect 
of quantum correlations out of equilibrium, or the role 
the steady states are few of the open questions. To ad- 
dress these questions, efficient probling tools are neces- 
sary [12]. 

Here we propose to probe dynamical correlations in 
strongly interacting ultracold gases by combining a quan- 
tum non-demolition matter-light interface with a quan- 
tum memory. The idea, stemming from quantum re- 



peaters, exploits the well known fact that light is a good 
carrier of information however difficult to store, while 
matter can be made to store coherently information for 
long times, serving thus as an efficient quantum mem- 
ory. Specially suited for such a goal is the QND Faraday 
interface which maps the spin-polarization of atoms into 
the polarization of photons. Such QND measurement 
has been recently implemented to spin squeeze states in 
cold atomic ensembles well below the projection noise 
limit [14]. It has also been presented as an efficient tool 
to probe quantum magnetism in ultracold atoms as pro- 
posed in [15-19] and experimentally demonstrated [20]. 
In our proposal, the QND interface is used both, to coher- 
ently map a relevant operator of the many-body system 
J at different times (t, t + r, . . . ) to the light with- 
out performing a measurement at each time, and also 
to store J into a quantum memory between two con- 
secutive matter-light interactions for a time r compa- 
rable to the time scale of the internal many-body dy- 
namics. A final single measurement which includes a 
readout of the memory yields dynamical spin correlations 
{^p\J{t)J{t + r)J{t + r') . . . Itp), for some given state tp. 
This quantum memory-assisted probing (QMAP) scheme 
is minimally destructive, does not depend on the linear 
response of the system to an external perturbation [13], 
and is equally suitable for systems out of equilibrium, 
i. e. when is not an eigenstate of the many-body Hamil- 
tonian. We show that the correlations obtained using 
QMAP are fundamentally different from those obtained 
by repeating a QND measurement of J^(t) at different 
instances of time and then performing a correlational 
analysis. This difference comes from the quantum in- 
terference effect present in the QMAP process, since a 
single measurement is performed after coherently storing 
the observable of interest at different times. 

Before explaining in detail our protocol, let us first 
review the basic elements of a Faraday QND interac- 
tion. When a strongly polarized light beam interacts 
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off resonantly with the internal spin degrees of freedom 
of an atomic system the polarization of the photons ro- 
tate, an effect known as the Faraday rotation (see the 
recent review [21] and reference therein). If the light 
beam is e.g. strongly polarized along the x-axis and 
propagates along the 2;-axis, it can be fully described by 
time-integrated canonical operators Xl = S2/ y^A^ph and 
Pl — ^s/^/Aph, where ^'2 (^3) is the Stokes operator 
corresponding to the difference in the number of photons 
in the ±45 (in the two circular) polarizations, and A^ph 
is the total number of photons in the beam. After the 
Faraday interaction, the integrated equations of motion 
result into [16, 21] 

XT' = Xi" - Kj{Q), (1) 

where X]^ {Xj^^) is the light quadrature before (after) 
the interaction The coupling strength can be expressed 
as /5: = ^/dr]A^ where d is the optical depth of the atomic 
sample, and r]A is the spontaneous emission probabil- 
ity induced by probing. An optimal QND interaction is 
achieved with tja in the range from ^ d~^^'^ for single 
color probing to the constant ^ 1/3 for two-color prob- 
ing [21]. For cold samples d > 100 [23] and hence /^^ ^ 10 
should be within reach. The observable 

corresponds to the total modulated magnetization of the 
atoms illuminated by the light beam. A^ denotes the to- 
tal number of atoms confined in an optical lattice, and 
the 2:-component of the atomic spin at site n. The 
modulation, given by the coefficients c^, reflects the spa- 
tial dependence of the light beam intensity; in a standing 
wave conflguration [16], Cn = 2 cos (/cna — a)^, where k 
is the wave number of the probing laser, a is the opti- 
cal lattice spacing, and a describes the shift between the 
probing standing wave and the optical lattice. 

The many body operator J is di QND observable [24] 
since it commutes with the Faraday Hamiltonian used 
for the measurement. As required to measure non-trivial 
dynamical correlations, J7 does not commute with the 
many-body Hamiltonian. The light-matter interaction 
time can be chosen in the /isec range [14], i.e., much 
shorter than the relevant timescale of the many-body 
system and be considered instantaneous, so that the 
QND character of the interaction is preserved. As shown 
in [19], J corresponds to a non-trivial observable of mag- 
netic systems. The aim of this work is to use the QND 
probing scheme to access to correlations functions as 
{tlj\J'{t)J'{0)\tlj) for any initial state tjj. We concentrate 
on two-time correlations to explain the method but ex- 
tensions to n-time correlations are straightforward. 

To better understand the role of the quantum mem- 
ory, let us first analyze the information acquired by just 
performing independent QND measurements of J at dif- 
ferent times (without storing them in a memory) and 




FIG. 1: (Color online) Schematic representation as a quantum 
circuit of the protocol to measure dynamical correlations. S 
(M) is the many body system (the quantum memory), and 
Li is the i-th light beam used in the protocol. The definition 
of the quantum gates is provided in the text, a) Scheme for 
measuring J' at two different instances of times, which after 
repetition give Fs{t). b) Scheme to measure FM(t) by using 
a quantum memory. 



making a statistical analysis of the results. Suppose our 
many-body system is initially prepared in some state 
At t = 0, a measurement oi J = "^-aiPi is per- 
formed, where Pi projects into the eigenspace with eigen- 
value equal to a^. Depending on the outcome a^, the 
state collapses to P^\^l^). In our set-up, this measure- 
ment corresponds to performing a QND Faraday inter- 
action between the system and the light, followed by a 
homodyne measurement of the relevant light quadrature. 
After the measurement, the many-body system evolves 
U{t) = exp[— zJ^MB^] under the many body Hamiltonian, 
Hub , during a time window of length t until a new mea- 
surement of J is performed yielding the outcome d j , as 
schematically sketched in Fig. la. The most general 
expression for the statistical mean resulting from both 
outcomes is given by 

Fsit)=J2a^{i'\P^J{t)P^\^)■ (3) 

i 

This quantity depends explicitly through expression (2) 
on the parameters k and a of the probing laser. If the 
state of interest is prepared in the ground state of the 
many-body system, i.e. = \Eq), the Fourier trans- 
form of Fs reads Cs{uj) = J dte'^^^Fsit) = Y^.^ ^ij6{uj - 
{Ei — Ej)), where {Ei} is the energy spectrum and the 
amplitudes (^ij can be trivially obtained. This quantity 
signals the frequencies corresponding to energy differ- 
ences between different eigenstates, providing partial in- 
formation about the energy spectrum. Although related, 
Fs{t) is not the two-time correlator of the observable 

To obtain two-time correlation functions we use optical 
quantum memories [25] based on atomic ensembles [21]. 
In the atomic memories, atoms, analogously to light, 
are strongly polarized along one direction, such that the 
spatially integrated spin components in the orthogonal 
plane, A'm(^) and PM(t), fulfill canonical commutation 
rules. The Faraday interaction between the light and the 
quantum memory yields the relation 

X<S}'=XZ+nwPT, (4) 
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where X]^ = XM{t = 0) and is the output value af- 
ter the interaction with the hght beam. The hght quadra- 
ture is Pl{z = 0). Note that both hght and mat- 
ter quadratures can be rotated according to X ^ X = 
cos(^)X + sm{(j))P and P P = cos((/))P - sin((/))X. 
The quadrature stored in the memory can be retrieved 
using a second hght beam and the corresponding Faraday 
interaction, such that X£^* = X)^ + i^rPm • Such type of 
memories have been experimentahy reahzed with atomic 
samples at room temperature in [26], where a storage 
time of the order of 4 ms was reported. Alternatively, 
memories with ultracold atoms in optical lattices have 
achieved longer storage times up to 240 ms [23]. These 
time storages are typically large enough to address the 
evolution of many-body ultracold gases. 

The protocol leading to spin dynamical correlations 
using a Faraday probing and a quantum memory 
— sketched in Fig. lb — , is straight foward and requires 
some rotations of the quadratures to obtain finally the 
correct time correlation. It is summarized as follows: 

(i) At t = 0, a Faraday interaction (F) between the 
many body system S and the first light pulse Li maps 
the (modulated) magnetization of the system J7 to the 
light quadrature according to X£^* = X)J^^ — /^iJ7(0) and 
Plf = P'j^^. A rotation of ^ = 7r/2 is applied to the 
output light quadrature. 

(ii) The rotated quadrature of this first light beam is 
mapped into the quantum memory M according to: 

xx?* = x» + KwPlT = ^iS - - «i J(0)). We 

denote this step by W "writing" in the scheme. 

(iii) The quadratures of the atomic memory are now 7r/2 



rotated; X^^ 
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(iv) The many-body system is let to evolve freely during 
an interval of time t such that J{t) = {t)JU{t). Due 
to the QND character of the Faraday interaction, is 
unaffected by the first interaction used to map its value 
to the memory. 

(v) At time t, a second light beam L2 is sent through 
both the system and the memory and reads (R) 

(vi) Finally, a balanced homodyne measurement of X£^*' 
is performed, denoted by D (detection) in Fig. lb. 

The variance of this observable is the key probe of the 
protocol, and is given by 



(5) 



where we have defined the noise of the signal as 



r]{t) = M^4[AJ{t)]^ 
i^ii^2f^Rf^w^ and A/" = (1 



^//^^[AJ(0)]^ being a^t 



)/2. Assuming as 



usually that the input laser field are in a coherent state 
and the memory is in coherent spin state, i.e. (X^") = 
(P^^) = and [AX^"]2 = [AP^"]^ = 1/2, the JV noise 
can be safely neglected, e.g. with K,i^hi2 ^ i^r^i^w > 1- 
By performing the two parts of the protocol indepen- 
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FIG. 2: (Color online) Fourier transform of Fm and Fs for 
Hamiltonian (Eq. (7)). The results are obtained with exact 
diagonalization for a chain of 12 spins, with a standing wave 
configuration such that k = 7v/{2a) and a = 0. Cm (filled 
red circles and dashed line) and Cs (empty blue squares and 
solid line) are plotted as a function of cj. Vertical dashed green 
lines show the spectrum of the Hamiltonian. The inset shows 
the strongly varying Fm (red), Eq. (6), and the almost flat 
Fs (blue), Eq. (3), as a function of time. The uncertainty in 
the QMAP signal due to experimental limitations is plotted 
as error bars (losses in the quantum memory lead to a 5% 
reduction of the signal), see the supplementary material for 
details. 



dently, i^l[AJ{t)]^ and 1^^ / kI[AJ (0)]^ can be measured 
and thus subtracted. After a tedious but elementary 
derivation, the signal obtained is 

FM{t) = mt), J{0)U) - 2{J{t)){J{0)). (6) 

which provides direct access to the symmetrized two- 
time dynamical correlation function of J7. Using eqn. (2) 



leads to FM{t) = Y^n^rn^nCm [Gmnit,{)) + Gnm(0,t)] /X, 

where the two point spin time correlation function is 
G™„(t,t') = {jm^nit')) - {fm{t)){jn(t')) and the CO- 
efficients Cn depend explicitly on a and k. 

If the state of interest is an eigenstate of the many- 
body Hamiltonian, the Fourier transform of Fm reads 
Cm{uj) (X EnU^i^ - {En - Eo)) + S{uj + {E^ - Eo))] 
which provides direct information of the energy spec- 
trum. Indeed, Cm{^) can be related to the symmetric 
spin dynamical structure factor [27], which is well known 
to give access to the energy spectrum and the dispersion 
relation of the system. In this case, r]{t) is constant and 
the noise becomes a peak at zero frequency and can be 
thus easily distinguished from the signal with features at 
finite frequencies. 

To illustrate the differences between Fm and Fs we 
study an array of coupled double well superlattices of 
2X spins: 



H 



N-l 



2n ' J2n+1 ~r ^2j2n+l ' j2n+2 



n=0 



(7) 
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where j = {j^^j'^^j^)^ j^^y^z ^^le usual spin-1/2 oper- 
ators, and ^1(2) is the couphng between even-odd (odd- 
even) spins. This Hamiltonian has been implemented 
with optical superlattices (see [6, 28] and references 
therein). In this case gi{2) — 4^i(2)/^' where t\{2) is 
inter (intra) double well hopping rate, and \J the on-site 
interaction energy, and ~ 10 ms which is of the same 
order of reported storing times (T) in quantum memo- 
ries, as cited before. The condition gT ^ 1 is indeed re- 
quired to resolve the dynamics of the many-body system, 
which demands either longer storage times or faster dy- 
namics. Remarkably, a recent experiment [8] simulating 
antiferromagnetism with tilted optical lattices provides a 
faster dynamical timescale given by only the tunneling 
rate, which relaxes the requirements on the memory. We 
address first the equilibrium case, where gi = g2 = g and 
the initial state is the many-body ground state |£^o)- In 
Fig. 2, the discrete Fourier transform of Fs and Fm us- 
ing exact diagonalization are plotted with the parameters 
given in the caption of the figure (the inset shows Fs and 
Fm)- The signal obtained with the statistical analysis Fs 
is very weak and contains many frequencies, therefore the 
Fourier transform Cs is almost fiat. With the use of the 
quantum memory, much less frequencies are present in 
the signal and Cm shows large peaks at some of the en- 
ergy levels, including the energy gap to the first excited 
state. The robustness of the QMAP method in the pres- 
ence of experimental limitations and memory loses can 
be inferred from the error bars in Fig. 2 as explained in 
the supplementary information. One could also study the 
non-equilibrium case by considering the system initially 
prepared in the ground state of the Hamiltonian Eq. (7) 
with gi = g and and ^2 = 0, that is, a product of singlets 
between odd-even spins. Then, at t = 0, the Hamiltonian 
is quenched to gi = g2 = g^ and the state would evolve 
accordingly. Our method could explore thermalization 
effects characterized by dynamical correlations [29], an 
interesting study which is beyond the scope of this work. 

In conclusion, we have shown that by combining a 
QND light-matter interaction and a quantum memory, 
a signal which exploits quantum interference provides di- 
rect access to dynamical spin-spin correlations. It is re- 
markable that in the presence of a many-body Hamilto- 
nian which does not commute with the QND interaction, 
even the storage of a single quadrature in the quantum 
memory yields results different from those obtained with 
a classical memory. The integration of quantum mem- 
ories in coherent spectroscopy techniques offers uncon- 
ventional possibilities to manipulate the dynamics of the 
many body systems by performing conditional feedback 
operations on the quantum memory [21, 30]. Particu- 
larly appealing are also prospects for detection of multi- 
time correlations or the study of non-equilibrium physics 
in strongly correlated systems. 
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